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Abstract. We investigate a class of dark energy models in which the equation of state undergoes a rapid transition and for which 
the Hubble SN la diagram is known to be poorly discriminant. Interestingly enough, we find that transitions at high redshift 
can lead to distortion in the correlation function of dark matter at lower redshift. We therefore use a combination of the SN la 
Hubble diagram, Cosmic Microwave Background data and power spectrum from the Sloan Digital Sky Survey Luminous Red 
Galaxies (SDSS LRG) to constrain the redshift of a possible transition. We find that the fundamental cosmological parameters 
are well constrained independently of the presence of a transition. Acceptable transitions from an equation of state close to 
w = to a value close to -1 are strongly rejected at redshifts much higher than those for which Large Scale Structure and 
SN la data are available: the transition redshift can be rejected up to a value as high as 10. We conclude that no preference for 
a transition appears from present-day data. 
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1. Introduction 

In 1998 and 1999 two groups collecting far supernovae data 
presented evidence for an acceleration of the expansion of the 
universe, and consequently for a non-zero cosmological con- 
stant, based on the Hubble diagram of the SN la supernovae 
dRiess et al. 19981 IPerlmutter et al. 1999b . Since these times 
published supernovae datasets have continuously increased and 
strengthened the case for an acceleration of the expansion of 
the universe. This has been corroborated by others data sets 
dFrieman et al. 20081 IBlanchard 20101 ). Although the cosmo- 
logical constant is entirely consistent with the detected accel- 
eration, the addition of this term to the theory of general rel- 
ativity has no known direct motivation. However, the cosmo- 
logical constant can also be interpreted as arising from the vac- 
uum contribution to the energy-momentum tensor. Indeed, the 
contributions of the quantum fluctuations of all the fields fill- 
ing the universe provide a non-zero density to the vacuum and 
a negative pressure acting exactly as a cosmological constant. 
However, such anticipated contribution from quantum fluctu- 
ations is many orders of magnitude larger than the presently 
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observed value. The low level of the energy scale of the cos- 
mological constant is therefore a mystery. Facing such a dif- 
ficulty in understanding this discrepancy, one can consider 
the existence of some unknown mechanism that cancels the 
contribution from vacuum energy, and look for other phys- 
ical components in the Universe that can produce an accel- 
eration of the expansion. Motivated by these considerations, 
|Ratra & Peebles (T 988) showed that the presence of a dynami- 
cal scalar field, not exactly at rest (otherwise, it would behave 
exactly as a cosmological constant) can actually generate an ac- 
celerated expansion. This scenario, known as the quintessence 
paradigm, has received a lot of attention in recent years af- 
ter the discovery of the acceleration. Other options have been 
proposed for its origin. An attractive idea is the possibility 
that acceleration appears as a non-linear contribution from in- 
homogeneities (Buchert 2008). However, no convincing argu- 
ments have been proposed to show that the actual contribution 
from these non-linearities can be brought to an observable level 
much above the naive expectation < h 2 >~ 10~ 10 , nor that 
there is reason why this contribution could have a significant 
apparent "negative pressure" to actually produce an acceler- 
ation. Another more radical option is to modify the equations 
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of the general relativity ( Cognola et al. 2006 ). However, adding 
new "exotic" components to the universe to obtain the acceler- 
ation is the most simple solution, and quintessence is natural in 
this context: in such models dark energy occurs from the late 
domination of some scalar field <p. The canonical Lagrangian 
of such a scalar field is given by 



-C = ^(dicf>) 2 - V(<f>) =X- V(cf>) , 



(1) 



with X being the kinetic energy of the field and V(4>) the poten- 
tial. The equation governing the evolution of this homogeneous 
field in an expanding universe is 



dV 
dtp 



(2) 



The stress-energy tensor has a form identical to that of an ideal 
fluid with pressure and density given by the two relations 



P = L = \(d4? ~ V((f>) 



p = 2 X£ x -£ = -{dS 2 + V(<P) , 



(3) 



(4) 



where x stands for — . For a field which is spatially homoge- 
dX 

neous the parameter w is defined from the equation of state 



p \j> 2 - V(<p) 



w = — = 



H 1 + v(4>) 



which remains greater than -1, while the sound velocity 
-C,x 



r 2 _ P* 



p x £x + 2X£ 



= 1 



(5) 



(6) 



AT 



is equal to the speed of light. However, more gen- 
eral Lagrangian for the scalar field have been advocated 
dArmendariz-Picon et al. 20001 lAvelino et al. 20091 1. In these 
models, known as ^-essence models, the Lagrangian is writ- 
ten in a very general form £ - p(X, (/>), allowing w to 
be smaller than -1 and eventually leading to an arbitrary 
value of Cj (lArmendariz-Picon et al. 200T1 ). Possible observa- 
tional consequences of a sound velocity c s different from the 
speed of light leads to effects which are below current limits 
(Hannestad 2005). Therefore, for simplicity, we assume below 
that the sound velocity remains equal to the speed of light and 
treat w(z) as a phenomenological function, which should be de- 
termined from observations. 

As the damping time scale in Eq.|2]is (3//) _1 , the kinetical 
energy of the field is not expected to change over a time scale 
much shorter than the Hubble time. Therefore, for smooth po- 
tentials, w(z) is expected to vary over a time scale typical of 
the Hubble time. For these slowly varying models, the CPL 
model (Chevalier & Polarski 200 1ILinder 2 003 ) is a useful ap- 
proximation that catches the essential feature of many models 
of varying dark energy: 



wcpl(z) = Wo + W\ 



1 + z 



(7) 



This model was used to define a figure of merit for a given set 
of astrophysical observations (Albrec ht et al. 200 6 ). However, 



rapid transitions in w at low redshift are not well described by 
such parametrization, and the correct dynamic of dark energy 
should be taken into account in cosmological parameter esti- 



mations ( Virey et al. 2005| l. It is therefore interesting to inves- 
tigate a wider class of parametrization. This is the purpose of 
the present work. In the second section we present the basics 
of the phenomenology of the quintessence scalar field and the 
parametrization used in this study. In the third section we ex- 
amine the constraints coming from three commonly used ob- 
servables, the type la supernovae Hubble diagram, the power 
spectrum of the cosmological microwave background (CMB) 
data published by the Wilkinson Microwave Anisotropy Probe 
(WMAP) collaboration and the power spectrum of Large Red 
Galaxies (LRG) from the Sloan Digital Sky Survey (SDSS). 

2. Dark energy parametrization 

For a pressureless regime like non-relativistic matter one has 

w = 0, for radiation w - -, while for a cosmological 

3 6 
constant w = -1. As we mentioned above, rapid transi- 
tions in w at non vanishing redshift are not well described 
by CPL parametrization ( Linden & Virey 2008) ). It is there- 
fore interesting to investigate a wider class of parametriza- 
tion. A convenient parametrization is (Lin der & Huterer 2 005 ; 
Douspis et al. 2008| > 



1 1 / 

w t (z) = -0+ + w-) - -(w+ - w_)tanhlFlogl 



(8) 



The equation of state parameter w, evolves from w,(z = +°°) = 
w+ to w t (z = -1) = w_, the transition occurs at redshift Zt or 
equivalently a t = 1/(1 + z,). The parameter Y figures the hard- 
ness of the transition: high values of F correspond to rapid tran- 
sitions, i.e. transitions with a shorter duration than the Hubble 
time atz,, while F ~ 1 corresponds to transitions occurring over 
a period similar to the Hubble time. The above parametric form 
may hide complexities : w < -1 (Caldwell 2002 ) models might 
be pathological dCarroll et al. 20031 . We try to keep our investi- 
gation as open as possible and do not add theoretical constraints 
on the kinds of models we investigate. For instance, modi- 
fied gravity theories might lead to an evolution in w(z) which 
crosses the phantom limit of w — -1 (Gannouji et al. 2006 



Bouhmadi-Lopez & Ferrera 2008 ). We therefore do not impose 



any restriction on the classes of models we investigate in terms 
of the values of w + and w_. 

With an appropriate choice of the parameter T, this 
parametrization is similar to the CPL parametrization (Eq. |7]l 
for z > 0: 



1 

wq = -(w+ + W-) 



1 

W\ = -(w+ - W-) 



(9) 



Requesting that the derivative of w, function at z = is the 
same as the derivative wcpl function (w',(z = 0) = w' CPL (z = 0)) 
imposes F = 1 and a, = 1 . However, slightly different values 
of T might lead to a w,{z) closer to w(z) over a wide redshift 
range. For instance, with F ~ 0.85, w(z) is very close to the CPL 
parametrization in the redshift range to 1 (see Fig. []]). The 
above parametrization provides therefore a nice generalization, 
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Fig. 1. Left: the evolution of the equation of state parameter of dark energy w(z) from Eq. [8] is illustrated on this plot for a 
transition redshift Zt = 1 : high values of F correspond to fast transitions, i.e. whose duration is much shorter than the Hubble 
time. Right: comparison between CPL parametrization and Eq.[8] 



comprising a close approximation of the CPL case, but also 
allowing the investigation of more rapid transitions. 

3. Constraints from present-day cosmological data 

Modifications of the equation of state of dark energy obvi- 
ously lead to modification of the expansion rate through the 
Friedmann-Lemaitre equation 



R = -4nG(p m +p x (l + 3w(z)) 



(10) 



This has direct consequences on geometrical tests. Another im- 
portant aspect comes from the fact that the growth rate of fluc- 
tuations in the linear regime is modified because the nature of 
the background is changed. For w+ ~ the contribution of 
dark energy can be non-negligible at high redshift and there- 
fore may leave specific imprints in the CMB and in the matter 
power spectrum. These imprints may not be reducible to a sin- 
gle geometrical factor (Ferram acho et al. 20091 1. 

3.1. Supernovae Hubble diagram 



In a preceding paper (Douspis et al. 2008) we presented some 
evidence that supernovae alone cannot reject a sharp transition 
between two values of w: w+ and w_ above a transition red- 
shift Zt(z) = 0.2 at 2 sigma using a compilation of supernovae 
recently published dDavis et al. 20 07 ). We chose for the two 
extreme values of w, -0.2 for w+ and -1 for w_. The choice of 
the second value was linked to the observation that the value 
of w now seems very close to -1, and the choice of w_ was 
discussed in detail in the article. There was also some surpris- 
ing but weak evidence that supernovae prefer a sharp transition 
at low redshift (Bassett et al. 2002b . The number of published 
SN la samples has increased in recent years. In 2008, Kowalski 
and coworkers (Kowal ski et al. 20 08 ) have applied a homoge- 
neous treatment to 13 samples of SNla published to reduce 



the possible systematics due to the analysis done under dif- 
ferent conditions and by different teams as much as possible. 
The tendency of SN la sample to prefer a transition has still be 
claimed to be present in this last homogeneous supernovae data 
set dShafieloo et al. 20091) . This provides us with motivation for 
more detailed investigations. 

3.2. Constraints from the CMB 

In order to use measurements of the CMB fluctuations to con- 
strain models with a rapid transition in dark energy we need 
to implement these models in a code which computes the 
CMB fluctuations angular spectrum. We therefore use a modi- 
fied version of CAMB in which our dark energy parametriza- 



tion has been implemented (see |Douspis et al. 2008 for de- 
tails on this implementation). As the value of the dark en- 
ergy density is analytically known in this model, this imple- 
mentation is straightforward, something which might not be 
the case if it has to be evaluated through numerical integra- 
tion (Corasanit i et al. 20041 1. as is generally the case. The sound 
speed was kept constant. We did not use the reduced quantity R, 
which is closely related to the angular distance to the CMB, and 
which is the most changed quantity when adding a cosmolog- 
ical constant (jBlanchard 19841 lEfstathiou & Bond 1999b . We 
used the 5-year release (INolta et al. 20091 1 and the ACBAR 
data (Reichar dt et al. 20 08 ) for additional stringent constraint 
on small scales. 



3.3. The effects of a transition model for dark energy 
in large scale structure data 

Before trying to constrain any fast transition model us- 
ing all data, it is interesting to examine how these mod- 
els affect the predicted matter distribution at low redshifts. 
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For this purpose, we used the correlation function data of 
Eisenstein et al. (2005) E05 hereafter, which provide a help- 
ful visualization of baryonic acoustic oscillations, since they 
appear as a well defined feature, the so-called acoustic peak in 
this correlation function corresponding to the oscillations in the 
power spectrum (BAO for baryonic acoustic oscillation) . 

We applied the procedure for determining the expected cor- 
relation function as described in Ferramacho et al. (2009), with 
a modified CAMB routine that comprises our varying dark 
energy parametrization (Douspis et al . 2008} . As discussed in 
Ferramacho et al. (2009), non-linear corrections and scale de- 
pendent bias play an important role in the final constraints 
obtained from LSS distribution and BAO. We have thus ap- 
plied the same prescriptions for these corrections, but taking 
some caution. For instance, we did not include the model- 
ing that accounts for the suppression of the baryonic acous- 
tic peak, since it uses the no-wiggle approximation from 
Eisenstei n & Hu 19981 which comprises fitting formulas de- 
rived for a standard ACDM cosmology. We checked however 
that correcting for this effect does not change the overall con- 
straints on w and the density parameters, except for Qi, which 
is expected since the amplitude of the acoustic peak depends 
mostly on the baryonic content of the universe rather than on 
any other quantity. In order to perform these corrections at 
low scales with less model dependence, it is also preferable to 
choose a more general correction model than the E05 approach. 
For this purpose, a Q-model (Col e et al. 20051 1 was used with b 
and Q as free parameters, which seems to be more appropriate 
for dark energy models with dynamical equation of state. 

In Fig. [2] we illustrate the effect of fast transitions in the 
correlation function at low redshift by considering three differ- 
ent values for w + and setting w_ = -1. In each plot, the effect 
of changing the transition epoch (a,) is shown. The remaining 
cosmological parameters were set to the best values obtained 
in Ferramacho et al. (2009) when considering a free but con- 
stant value for w. The bias parameter was marginalized over in 
all the plots for a better comparison of the effects in the over- 
all shape. Matter fluctuation amplitude, measured by <x 8 , might 
therefore be different for the different models. 

Transitions between w + = -0.2 and w_ = - 1 do not change 
the overall shape of the correlation function, even if we con- 
sider transitions at relatively low redshifts. In this case, the 
predicted matter distribution with the adopted model for dark 
energy cannot be distinguished from the same distribution ob- 
tained with a standard ACDM model. However, when we let 
w + be closer to 0, some interesting effects appear in the pre- 
dicted correlation function, noticeably at intermediate and large 
scales. We observe then a relative increase of the power at these 
scales, which boosts the amplitude of the baryonic peak (to a 
small degree when w+ = -0.1 and much larger when w + - 0) 
and shifts its position, and furthermore the global shape of the 
correlation function is modified. These effects show that the 
matter distribution is sensitive to fast and strong variations in 
the dark energy equation of state and provides in principle a 
good probe to test such models. Also, in the view of these re- 
sults changes in the correlation function for these models can- 
not be fully reproduced with a single geometrical factor, as for 
instance the A parameter from E05. So, the full shape of the 
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Fig. 2. Expected shape for the correlation function in units of 
%r 2 , assuming rapid (step-like) transitions with r = 5 in the 
dark energy equation of state at different redshifts, with w_ = 
— 1. Upper panel: w+ — —0.2. Middle panel: vv+ = —0.1. Lower 
panel: w+ = 0. Plotted curves correspond to different values for 
the transition redshift: z, = 0.5 (dashed-dot), z t — 1 (dashed), 
Zt — 2 (dot) and Zt — +°° (straight). 



correlation function (or of the matter power spectrum) must be 
taken into account when studying models that have an effect on 
both the correlation function form and the position of the bary- 
onic peak due to other effects than angular diameter distance 
changes. The precision of the SDSS LRG measurements on the 
correlation function is not high, but as an example the model 
with w + = -0. 1 and a transition redshift fixed at z r = 0.5 has 
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Fig. 3. 2D contours at 68 % and 95% confidence level, for the 
parameters wq = l/2(w + + w_) and w\ — l/2(w + - w_), set- 
ting tt = 1. The blue contours correspond to T = 1 and the 
red contours to F = 0.85. The green countour shows the same 
confidence regions using the standard CPL parameterization. 

a x 1 reduced of 2 in comparison with the concordance model, 
suggesting that a transition is preferred at 1 <x using only the 
correlation function data. 

That transitions in w at redshifts as high as 1 could leave 
an imprint on the correlation function is an interesting issue 
requiring a further investigation. Below, we look into the im- 
pact on these models of the joint analysis of data from galaxy 
distribution, CMB anisotropics and SN la. 



a constant value (no evolution). The constraints on nearly all 
parameters remain essentially identical. Of course, the main 
exception regards parameters that govern the equation of state 
itself, which have a bigger uncertainty compared to a constant 
w. Considering the parametrization of a smooth transition with 
either r = 1 or F = 0.85 results in some small changes in 
the preferred values and the uncertainties, which seem to be 
smaller when T = 0.85. The 2D marginalized contours at 68% 
and 95 % confidence level are presented in Fig. [3] Again, the 
slightly different values of the F parameter have an effect on 
the derived constraints on wo - W\, although it is rather small. 
Finally, we notice that using the parametrization given by Eq.[8] 
with T = 0.85 provides very similar constraints to a pure CPL 
parametrization. 

The Dark Energy Task Force report ( Alb recht et al. 2 006 ) 
introduced a numerical quantity to determine the capacity of 
future surveys in constraining dark energy models. This numer- 
ical quantity is referred to as the figure of merit (f.o.m.) for a 
given experiment or combination of experiments and is defined 
as the reciprocal of the area of the marginalized 2D region con- 
tour at 95 % confidence in the wq— w\ parameter space. In terms 
of the f.o.m. for the present day data, the obtained values are 
summarized in the following table 



Model CPL r = 1 T = 0.85 
f.o.m. 2XB L54 2.11 



3.4. Hubble time transitions - CPL approximation 

Let us start our analysis by taking the case of transitions oc- 
curring at a rate typical of the Hubble time. As seen above, 
with F ~ 1, the parametrization provided by Eq. [8] becomes 
almost equivalent to the CPL parametrization. It is interest- 
ing to compare the effect of considering r = 1, which fulfills 
the condition w t (z = 0) = W\, as requested when deriving the 
CPL parametrization, and a slightly lower value for F like 0.85, 
which provides a better agreement for the overall shape of the 
w(z) curve in both parametrizations. 

Both cases were evaluated together with the CPL 
parametrization (Eq. [7]). We performed a MCMC run 
by introducing the modified CAMB routines and the 
new parametrization for w(z) in the COSMOMC code. 
The total parameter space directly constrained was then 
p={Q/,/z 2 , Q. c h 2 , 6, t,A s , n s , wo,wi}. We used the matter power 



spectrum from Tegmarketal. 2006 to account for the SDSS 
LRG data. The modeling for non linear effects and scale de- 
pendent bias was done with a Q-model dCole et al. 20051 1 as 
explained above. Since we are interested in models with w{z) < 
- 1 in order to determine a valid figure of merit, dark energy 
perturbations were "turned off" when solving for CMB and 
matter power spectrum. 

The obtained results on all parameters are summarized in 
Table [TJ All these constraints on individual parameters were 
obtained by marginalization over other parameters. For com- 
parison reasons, the last column shows the results where w is 



These values can thus be used to compare the improvement 
of future experiments in constraining dark energy evolution on 
Hubble timescales. For more rapid transitions, one has to con- 
sider higher values of T, which is done in the next section. 



3.5. Joint constraints 

In this section we examine what constraints can be derived 
when a rapid transition in the equation of state of dark energy 
(r = 5) is considered, using the full data sets presented above. 
Only models with w(z) > -1 were constrained by including 
DE perturbations when solving the total perturbations equa- 
tions with CAMB. For each case we ran one long chain, us- 
ing the Rafetery-Lewis diagnostic to check for convergence. As 
for the LRG SDSS data, we used the data on the matter power 
spectrum as implemented in COSMOMC with a Q-model to 
correct for non-linear effects. Finally, we restricted ourselves 
to a flat Universe (Q^ = 0) for simplicity reasons. 

As seen above, for transitions to w_ = - 1 , the matter corre- 
lation function shows a different sensitivity to the initial value 
w+. We considered then three distinct transitions by setting 
W- = — 1 and three low values for w+ (0,-0.1 and -0.2). In 
each case the objective is to get constraints on the transition 
epoch (a,), the value a, = corresponding to the standard con- 
cordance model. As w+ approaches the contribution of dark 
energy at high redshifts remains non vanishing, allowing us to 
investigate its possible influence at earlier epochs. 
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3.5.1 . • W- — — 1,M>+ = -0.2 

Four different data combinations were used to constrain 
this model, CMB, CMB+SN la, CMB+P(k) and finally 
CMB+P(k)+ SN la. The results of the mean posterior values 
and 68% confidence intervals on each parameter are summa- 
rized in Table [2] The posterior distribution in some parameter 
is clearly non-Gaussian, so in some cases the lcr error bars are 
quoted using asymmetrical limits. Fig.|4]shows the contours in 
the space of Q.q - a t , on which the effect of the different data 
combinations on these two parameters can be appreciated. This 
figure shows that the combination of the WMAP5 data on the 
CMB with supernovae already constrains the possible transi- 
tion to be at a, < 0.6 at the two sigma level (95% confidence 
region), or in redshift space z, > 0.7. When the matter power 
spectrum data is added, some small reduction on the allowed 
region in the transition parameter a, is observed. Surprisingly, 
the combination CMB+P(k) without supernovae leads to a bet- 
ter constraint on the transition epoch. This result is not intuitive, 
since as demonstrated in the previous section, one obtains only 
negligible differences in the matter correlation function shape 
by changing the transition epoch. The cause for this reduction 
of the allowed transition redshift is actually due to the shift in 
the preferred Q„„ as the power spectrum data tend to prefer 
lower values of Q m . 

Another interesting factor of this dark energy modeling 
is the effect that transitions at intermediate redshifts have in 
the allowed matter fluctuation amplitude - the cr% parameter. 
Figure [5]presents a peculiar degeneracy between these param- 
eters. This figure shows that for transitions occurring between 
a, * 0.2 (z t = 4) and a, « 0.6 (z, = 0.6), the allowed mat- 
ter amplitude cr 8 decreases linearly, meaning that transitions 
which occur at lower redshifts necessarily yield a lower value 
of erg. This is because a value of w close to has the effect 
of inhibiting the growth of matter perturbations, as shown in 



Douspis et al. 2003| So, in the presence of a transition from a 



value of w close to 0, the overall amplitude of fluctuations (for 
a given CMB amplitude) is smaller at low redshifts. In princi- 
ple the amplitude of matter fluctuations can be estimated from 
local cluster abundances. However, in practice the lack of reli- 
able calibration of the mass temperature limits the use of this 



constraint (Pier paoli et al. 200T]|Blanchard & Douspis 2005[ ). 



3.5.2. w- = -l,w+ = -0.1 

We have seen above that for a value of w+ closer to impor- 
tant changes in the predicted matter correlation function arise, 
which can improve the quality of the fit of the SDSS LRG data. 
We now check whether these models, which comprise transi- 
tions at low redshifts, may accommodate other observational 
data. So, we determine the same combined constraints as done 
previously, but setting w+ = —0.1. The obtained results for the 
other data combinations are summarized in Table [3] and in Fig. 
|6l which shows the marginalized region contours in the Q.q - a, 
parameter space. For these models the allowed transition red- 
shift space is reduced: models with Qg 0.75 (Q m 0.25) 
and a, > 0.5 (z, < 1) are excluded at the two sigma level. This 
means that for these models the CMB data impose a stronger 
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Fig. 4. Confidence regions at 68 % and 95 % level for the tran- 
sition epoch a, and the density of the quintessence field Og, 
considering w+ - -0.2 and w_ = -1 and using different 
combinations of data:CMB {light blue), CMB+ SN la {red), 
CMB+P{k){green) and CMB+SN Ia+P(k)(W«e). 
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Fig. 5. Confidence regions at 68 % and 95 % level for the tran- 
sition epoch a, and the amplitude of matter perturbations erg, 
considering vv+ = -0.2 and W- = — 1 and using different com- 
binations of data: CMB+ SN la {red), CMB+P{k)(green) and 
CMB+SN Ia+P(k)(Wwe). 

constraint, and the best fit model in middle panel of Fig.|2]with 
Zt = 0.5 becomes ruled out by the CMB. 



3.5.3. w- = -l,w + - 

Finally we consider the case where the dark energy has a van- 
ishing equation of state parameter in the early Universe. This 
case is phenomenologically interesting as dark energy then 
evolves as pressureless matter, although it cannot be accounted 
for in standard quintessence models like the inverse power law 
or SUGRA potentials ( |Douspis et al. 2008| >. In Tableland Fig. 
Qwe present the same information as above for the constraints 
on this model. The result is somehow expected from what we 
saw when we passed from w+ - 0.2 to w+ = -0. 1 . While the 
constraints on Q. m (Qg) remain almost unchanged, the possi- 
ble transition epoch is further pushed towards higher redshifts, 
with a value Zt > 2 at 95% confidence level. Again, this con- 
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Fig. 6. Confidence regions at 68 % and 95 % level for the tran- 
sition epoch a, and the density of the quintessence field Q.q, 
considering w + = -0.1 and W- = — 1 and using different com- 
binations of data: CMB+SN la (red), CMB+P(k)(green) and 
CMB+SN la+P(k)(blue). 

straint comes mainly from the combination of the CMB obser- 
vations with the other probes. 

From the three considered cases, the models with late tran- 
sitions Zt < 1 seem to be excluded since they cause the angular 
distance to the CMB to change significantly which disagrees 
with the position of the acoustic peak in the temperature power 
spectrum. This finding does not allow for lower values of the 
transition redshift when combined with the strong constraint 
on Q„, from the SDSS LRG distribution at low redshifts and 
the SN la data. 



w 1 at low redshifts and in the future. This assumption was 

based on the theoretical prior that a final w = -1 is easier to 
understand in existing quintessence models. Nevertheless, one 
can relax this assumption on a purely phenomenological basis, 
and check what happens to the constraints derived above when 
W- is left as a free parameter. We have performed this analysis, 
taking the same three values for w+ as before and using the full 
combination of data CMB+P(k)+SN la. 

The results on the three chains are summarized in Table [5] 
The main constraints on the transition parameters are summa- 
rized in Fig. [8] and [9] which show the effect of leaving w- as 
a free parameter in the transition epoch defined by a,. The ef- 
fect is drastic, pushing the possible transition to much higher 
redshifts while at the same time the preferred value for w- be- 
comes w- - -0.956, with w_ = — 1 just compatible at one 
sigma. That the constraint put on a, for the different values of 
vv+ becomes much more severe (z t > 4 at two sigma, in the 
best case) may seem surprising since the obtained value is not 
so far from the fiducial w = -1 considered above. As a check, 
we reevaluated our constraints setting w_ = -0.95 and found 
identical constraints to those marginalized over w. 

It seems that by letting the present day value for the equa- 
tion of state be a little higher than -1, current data tend to 
strongly prefer the case of no transition, rejecting any accept- 
able transition to higher redshifts. This agrees with the con- 
straints obtained when we considered for a constant value of w 
previously, which is a good cross-check of this analysis. That 
the value w = -0.95 is systematically preferred is somehow 
interesting, although with the current constraint precision, the 
ACDM is obviously compatible with all derived constraints. 
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Fig. 7. Confidence regions at 68 % and 95 % level for the tran- 
sition epoch a, and the density of the quintessence field D.q, 
considering w+ — 0. and W- = — 1 and using different com- 
binations of data: CMB+SN la (red), CMB+P(k)(green) and 
CMB+ SN la+P(k)(blue). 



0.8 



0.6 



O 0.4 - 



0.2 



0.0 



0.5 



0.6 




0.9 



1.0 



Fig. 8. Confidence regions at 68 % and 95 % level for the tran- 
sition epoch a, and the density of the quintessence field Og, 
w_ being an additional parameter on which the likelihood is 
marginalized and using different values for w+: (green), -0.1 
(red) and -0.2 (blue) using CMB+SN Ia+P(k). 



3.5.4. Relaxing the assumption of w_ = -1 

We have seen the effects of considering transitions that have 
a final value W- compatible with a cosmological constant, i.e. 



4. Discussion and conclusions 

We have investigated quintessence models which undergo a 
transition in the equation of state parameter w of their dark en- 
ergy component in the light of latest cosmological data. We 
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Fig. 9. Confidence regions at 68 % and 95 % level for the transi- 
tion epoch a, and w_, marginalizing on Qg with different values 
for w+: (green), -0.1 (red) and -0.2 (blue). 



concentrated on a class of models in which the duration of 
the transition is much shorter than the Hubble time, for which 
we have previously shown that the supernovae Hubble diagram 
alone provides, somewhat surprisingly, poor constraints. A ma- 
jor difference to previous studies of varying dark energy models 
is that we fully took into account the data on the matter power 
spectrum rather than using a reduced geometrical factor A and 
used a class of models in which the dark energy component at 
high redshift does not cross the matter regime, i.e. w remains 
below 0. We concentrated on models in which the high red- 
shift equation of state parameter approaches zero (specifically 
w = -0.2, w = —0.1, w = 0.), without marginalizing because 
of the sensitivity to w when it approaches zero. The first notice- 
able result is the stability of cosmological constraints obtained 
for the fundamental parameters: all parameters that our five ta- 
bles have in common are nearly unchanged in comparison with 
those obtained for a constant w (Ferram acho et al. 20091 ). In 
this respect we noticed that a value of w = -0.95 + 0.05 re- 
mains preferred at the one sigma level. The amplitude of mat- 
ter fluctuations reveals some sensitivity to dark energy dynam- 
ics, as the erg parmeter is lower in the presence of dark energy 
with an equation of state parameter w > — 1, a result consistent 



with previous findings for constant w (Do uspis et al. 2008) . A 
second conclusion is that a transition is generally severely re- 
jected when combinations of data are used, although a transi- 
tion would be (weakly) preferred by a single set of data like 
supernovae or large scale structure. We notice however that the 
redshift at which the possible transition is rejected delicately 
depends on the assumed hypothesis: with a present-day w set to 
1 and w+ set to -0.2 we could formally reject at 2<x a tran- 
sition below z ~ 1 ■ 1 , while when we marginalize on w- and set 
w+ to 0, the transition is allowed only at redshifts greater than 
~ 10. Very recently, new data and analysis have appeared in the 
literature on SN la dFreedman et al. 20091 iKessler et aT~2 009 ; 



Lampeitl et al. 2009 1 and power spectrum ( Perci val et al. 2 009 ; 
Reid et al. 20091 1. It is therefore legitimate to wonder whether 



reveals that our results are not sensitive to different data com- 
binations: using the three data sets (CMB, P(k), SN la) or us- 
ing only two (either CMB and P(k) or CMB and SN la) leads 
to nearly identical constraints, differences in preferred param- 
eters being at most 1.5 <x, less than the worse systematic dif- 
ference discussed below. Since preferred cosmological param- 
eters from most recent analysis are entirely consistent with our 
previous analysis (Ferramac ho et al. 20091) . it is reasonable to 
assume that using the new data would make virtually no differ- 
ence. Our general conclusion is, maybe not surprisingly, that 
the concordance model is doing well and that possible varia- 
tion in dark energy is unlikely. 

A final word of caution though: In the present analy- 
sis we used the constraint from the power spectrum of LRG 



(Teg mark et al. 2006| ). However, Ferramacho et al. (2009) no- 
ticed that the use of the correlation function could lead to 
preferred values which differ up to 2cr from the ones ob- 
tained with the power spectrum. Therefore, our 2<x ranges 
can be regarded as indicative of possible systematic effects. 
We have no certainty however that this systematic difference 
could not be larger for varying dark energy models. The differ- 
ence in the clustering properties between red and blue galax- 
ies is also a subject that calls for caution in the detailed inter- 
pretation of data dSanchez & Cole 20081 IPercival et al. 20071 ). 
even if recent analysis tends to weaken earlier discrepancies 
(Perc ival et al. 20091 . Improvements of constraints on dark en- 
ergy from future experiments would have to control in an ex- 
tremely accurate way their systematics to infer robust con- 
straints on the variation of dark energy. 
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Table 1. Table resuming the obtained constraints for a pure CPL model (1), parameterization|8]with r = 1 (2), T = 0.85 (3) and 
assuming a constant w (4). 



Param. 


1 


2 




3 




Vanilla + w 


n b h 2 w 


0.0222 ± 0.0006 


0.0228 ± 


0.0006 


0.0228 =i 


: 0.006 


0.0228 ± 0.0006 


a c h 2 [»i 


0.107 + 0.005 


0.107 ± 


0.005 


0.107 ± 


0.005 


0.109 ±0.005 


Q\c\ 


0.104 ±0.003 


1.041 ± 


0.002 


1.041 ± 


0.002 


1.042 ±0.003 


T M 


0.856 ± 0.017 


0.088 ± 


0.017 


0.088 ± 


0.016 


0.087 ± 0.017 


n s w 


0.958 ±0.014 


0.969 ± 


0.008 


0.968 ± 


0.008 


0.967 ±0.014 


log(10 l0 A, s )W 


3.05 ± 0.04 


3.06 ± 


0.04 


3.06 ± 


0.04 


3.06 ± 0.04 




-1.027 ±0.123 


-0.994 d 


= 0.141 


-1.005 ±0.119 


-0.965 ± 0.056 


Wl [sl 


0.278 ± 0.408 


0.036 ± 


0.452 


0.186 ± 


0.404 





Q e (£2 A )[« 


0.735 ±0.017 


0.747 ± 


0.018 


0.738 ± 


0.018 


0.739 ±0.014 


Age [/1 


13.9 ±0.1 


13.7 ±0.1 


13.8 d 


:0.1 


13.7 ±0.1 


Q. W 


0.265 ± 0.017 


0.253 ± 


0.018 


0.262 ± 


0.018 


0.261 ± 0.020 




0.766 ± 0.037 


0.780 ± 


0.037 


0.776 ± 


0.037 


0.816 ±0.014 




11.0 ± 1.5 


11.1 d 


: 1.4 


11.1 d 


: 1.4 


11.0 ±1.5 




0.698 ±0.018 


0.718 ± 


0.020 


0.706 ± 


0.017 


0.713 ±0.015 



M Normalized baryon matter density times h 2 , m normalized CDM density times h 2 , w ratio of sound horizon to angular diameter distance, 
M optical depth to reionization, w primordial spectral index, W primordial scalar power at k = 0.05Mpc~', [;1 parameters describing the 
equation of state of dark energy, M dark energy density, ra age of the universe in Gyr, ^ normalized matter density, [ * ] matter density 
fluctuation amplitude, ra reionization redshift, [m] reduced Hubble parameter. 

Table 2. Table resuming the obtained constraints different combinations of data and setting w+ = -0.2 and w- = — 1. 



Param. 


CMB + SN la 


CMB + P(k) 


CMB + P(k) + SN la 




0.0225 ± 0.0007 


0.0224 ± 0.0006 


0.0224 ± 0.0006 


Cl c h 2 


0.107 ±0.007 


0.104 ±0.005 


0.105 ±0.005 


e 


1.040 ±0.003 


1.040 ±0.003 


1.039 ±0.003 


T 


0.087 ± 0.019 


0.088 ± 0.017 


0.085 ±0.016 


n s 


0.969 ±0.018 


0.967 ±0.016 


0.968 ± 0.017 


log(W w A s ) 


3.05 ± 0.04 


3.04 ± 0.04 


3.04 ± 0.04 


a, 


031+0.09 

-0.31 


25 +0 08 


rv 9Q+0.09 
"•^^-0.29 


n e (n A ) 


0.718 ±0.023 


0.743 ± 0.020 


0.733 ±0.017 


Age 


14 ± 0.2 


13.9 ±0.2 


14.0 ± 0.2 


o. m 


0.282 ± 0.023 


0.257 ± 0.020 


0.267 ± 0.017 


0-8 


0.705 ± 0.086 


0.717 ±0.063 


0.707 ± 0.075 


Zte 


11.1 ± 1.6 


11.1 ± 1.5 


11.0 ± 1.4 


h 


0.678 ± 0.028 


0.704 ± 0.025 


0.691 ±0.024 



Table 3. Table resuming the obtained constraints using different combinations of data and setting w + = -0. 1 and vv_ = - 1 . 



Param. CMB + SNIa CMB + P(k) CMB + P(k) + SN la 



a b H 2 


0.0224 ± 0.0007 


0.0224 ± 


0.0006 


0.0224 ± 0.0007 


n c h 2 


0.109 ±0.006 


0.105 ± 


0.005 


0.106 ±0.005 


e 


1.040 ±0.003 


1.040 ± 


0.003 


1.039 ± 0.003 


T 


0.089 ± 0.017 


0.089 ± 


0.017 


0.088 ±0.018 


n s 


0.970 ± 0.022 


0.971 ± 


0.019 


0.972 ± 0.020 


log{W°A s ) 


3.06 ± 0.04 


3.05 ± 


0.04 


3.04 ± 0.04 


a, 


25 +0 08 


23 +0 07 


"■ zu --0.26 




0.719 ±0.023 


0.743 ± 


0.020 


0.733 ± 0.017 


Age 


14 ± 0.2 


13.9 d 


:0.2 


14.0 ± 0.2 




0.281 ±0.023 


0.257 ± 


0.020 


0.267 ± 0.017 


0"8 


0.719 ±0.086 


0.708 ± 0.074 


0.696 ± 0.083 


Zre 


11.4 ± 1.5 


11.1 d 


: 1.5 


11.2 ± 1.6 


h 


0.686 ± 0.026 


0.704 ± 


0.024 


0.693 ± 0.023 
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Param. 


CMB + SN la 


CMB + P(k) 


CMB + P(k) + SN la 


Q. b h 2 


0.0223 ± 0.0006 


0.0223 ± 0.0006 


0.0223 ± 0.0006 


£l c h 2 


0.112 ±0.005 


0.108 ±0.004 


0.110 ±0.004 


e 


1.038 ±0.003 


1.038 ± 0.003 


1.038 ± 0.003 


T 


0.086 ± 0.016 


0.088 ± 0.017 


0.087 ± 0.017 


n s 


0.965 ±0.016 


0.968 ± 0.016 


0.965 ± 0.015 


log(10 10 A s ) 


3.06 ± 0.04 


3.05 ± 0.04 


3.05 ± 0.04 


a, 


16 +005 

-0.16 


15 +005 

-0.15 


16 +005 

-0.16 


n e (n A ) 


0.723 ±0.021 


0.746 ± 0.020 


0.737 ± 0.016 


Age 


13.8 ±0.1 


13.8 ±0.1 


14.0 ±0.1 


Q. m 


0.277 ±0.021 


0.254 ± 0.020 


0.263 ± 0.016 


OS 


0.759 ± 0.052 


0.740 ± 0.050 


0.746 ± 0.048 




11.1 ± 1.4 


11.1 ± 1.5 


11.1 ± 1.5 


h 


0.698 ±0.018 


0.717 ±0.019 


0.707 ± 0.015 



Table 5. Table resuming the obtained constraints using the full combination of data (CMB+P(k)+SN la) and considering w- as 
a free parameter. 



Param. 


w + = 


w + = -0.1 


w + = -0.2 


a b h 2 


0.0222 ± 0.0006 


0.0223 ± 0.0006 


0.0223 ± 0.0006 


0. c h 2 


0.106 ±0.005 


0.106 ±0.005 


0.107 ±0.005 


e 


1.039 ± 0.003 


1.039 ±0.003 


1.042 ± 0.003 


T 


0.085 ±0.017 


0.085 ± 0.016 


0.085 ± 0.017 


n s 


0.969 ± 0.008 


0.960 ±0.013 


0.960 ± 0.014 


log(l0 l0 A s ) 


3.05 ± 0.04 


3.05 ± 0.04 


3.05 ± 0.04 


a, 

W- 


-0.954 ±0.051 


q 7 j +0.018 

"• u/1 -0.071 

-0.951 ±0.058 


100 +0026 
-0.965 ± 0.056 


a e (n A ) 


0.742 ± 0.015 


0.742 ±0.015 


0.740 ± 0.014 


Age 


13.8 ±0.1 


13.8 ±0.1 


13.8 ±0.1 




0.258 ± 0.015 


0.258 ±0.015 


0.260 ± 0.015 




0.763 ± 0.036 


0.761 ± 0.038 


0.764 ± 0.037 


Z re 


11.0 ± 1.5 


11.0 ± 1.4 


10.9 ± 1.5 


h 


0.705 ± 0.016 


0.705 ±0.016 


0.704 ± 0.016 



